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ABSTRACT 

We present a new method to test the cosmological model, and to estimate the cosmological param¬ 
eters, based on the non-linear relation between ultraviolet and X-ray luminosity of quasars. We built 
a data set of 1,138 quasars by merging several literature samples with X-ray measurements at 2 keV 
and SDSS photometry, which was used to estimate the extinction-corrected 2500 A flux. We obtained 
three main results: (1) we checked the non-linear relation between X-ray and UV luminosities in 
small redshift bins up to z ~ 6, confirming that it holds at all redshifts with the same slope; (2) we 
built a Hubble diagram for quasars up to z ~ 6, which is well matched to that of supernovae in the 
common z = 0 — 1.4 redshift interval, and extends the test of the cosmological model up to z ^ 6; (3) 
we showed that this non-linear relation is a powerful tool to estimate cosmological parameters. With 
present data, assuming a ACDM model, we obtain r2M=0.22to;j8 and HA=0.92to;3o (Hm=0.28±0.04 
and Ha=0.73±0.08 from a joint quasar-SNe fit). Much more precise measurements will be achieved 
with future surveys. A few thousands SDSS quasars already have serendipitous X-ray observations 
with Chandra or XMM-Newton, and at least 100,000 quasars with UV and X-ray data will be avail¬ 
able from the eROSITA all-sky survey in a few years. Euclid^ LSST, and Athena surveys will further 
increase the sample size to at least several hundred thousands. Our simulations show that these 
samples will provide tight constraints on the cosmological parameters, and will allow to test possible 
deviations from the standard model with higher precisions than available today. 

Subject headings: Cosmology — Quasars 


1. INTRODUCTION 

An optimal “standard candle” for cosmological stud¬ 
ies is defined by two fundamental properties: it has a 
standard (or standard-izable) luminosity, and it is easy 
to observe in a wide redshift range. Quasars are the 
best class of astrophysical sources concerning the latter 
property, but completely lack the former: their observed 
emission spans several orders of magnitudes in luminos¬ 
ity, but their spectral energy distribution (SED) show 
no, or little, significant evolution with luminosity. As a 
consequence, quasars at first glance are far from being 
useful tools to set the cosmic distance scale as a function 
of redshift (i.e. the ’’Hubble Diagram”). 

However, even weak correlations between spectral fea¬ 
tures and luminosity, with large dispersions and obser¬ 
vational biases, can in principle be useful for cosmolog¬ 
ical measurements, provided that the quasar sample is 
large enough. This consideration prompted several at¬ 
tempts to derive cosmological parameters from quasar 
observations. Examples are the an ti-correlation between 
UV emission lin es and luminosity (Baldwin 1977 Osmer 


& Shields 19991, the lumino sity-mass relation m super- 
Eddington accreting quasars (Wang et al.|2014 1, the rela¬ 
tion betwee n luminosity and X-ray variability (L a Franca 
et al. 1120141 ) , the radius - luminosity relationship ([Watson 
et al. |2Ullf Melia||2014l |Kilerci Eser et al.||2015|) . 

In this paper we will explore the possibility of build¬ 
ing a Hubble diagram for quasars using the well known 
non-linear relation between the UV and X-ray luminosi¬ 
ties. In recent years, several quasar samples have been 
observed at both optical/UV and X-ray wavelengths. We 
will show that if we collect all the available data in the 
literature, we can build a large enough sample to signifi¬ 


cantly constrain the cosmological parameters, and to test 
the cosmological model over the whole redshift range out 
to z = 0 — 6.6 (age=0.8 Gyr). This analysis is achiev¬ 
able thanks to the increase of the number of quasars with 
observed UV and X-ray emission (provided by recent op¬ 
tical and X-ray surveys) by about an order of magnitude. 

A non-linear relation between the UV and X-ray lumi- 


X-ray surveys (Tananbaum et al. 111979 

Zamorani et al. 

1981 Avni &€ lananbaum|1986 

), and has been conhrmed 


with the main X-ray observatories over a redshift range 
from 0 to 6.5 and about five decades in UV luminosity. 
The largest samples recently analyzed include (1) a com¬ 
pilation of 367 quasars from different optical surveys and 
obser ved by ROSAT, XM M-Newton and Chand ra (33 3 
from Steffen et al. 2006 and 34 from Just et al. 2007); 
(2) a sample of 350 sources obtained from the Sloan Digi- 
tal Sky Surve y second quasar cata log with XMM-Newton 
observations (|Young et al.||2010|), (3) a sample of 545 
object from the GUSiVKJS-XMM-Aewton survey (Lusso 


et al. 2010); (4) a sample of 200 quasars with UV and 


X-ray observations from the Sw ift observatory (Grupe 
et al. 2010| Vagnetti et al.|[2M0 ). 


In all these works the Lx-Luv relation is parametrized 
as a linear dependence between the logarithm of the 


monochromatic luminosity at 2500 A (Luv) and the aox 
parameter, defined as the slope of a power law connect¬ 
ing the monochromatic luminosity at 2 keV (Lx), and 
Luv- aox=0.384x log(Lx/Luv)- Luminosities are de¬ 
rived from fluxes through a luminosity distance calcu¬ 
lated adopting a standard ACDM model with the best 
estimates of the cosmological parameters Dm and Da at 
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the time of the publications. When expressed as a rela¬ 
tion between X-ray and UV luminosities, the aox-Luy 
relation becomes: log Lx = /? + ylogLuv, with all the 
works cited above providing consistent values of the free 
parameters: /3 ^9 and 7 ~0.6. The observed dispersion 
is S ~0.35-0.40. 

The potential use of this relation as a cosmological 
probe is obvious: if we assume no redshift evolution of 
the relation, the observed X-ray flux is a function of the 
observed UV flux, the redshift, and the parameters of the 
adopted cosmological model. The relation can be then 
fitted to a set of UV and X-ray observations of quasars 
in order to estimate the cosmological parameters. 

However, none of the samples published so far has the 
size and/or homogeneity in the observations to provide 
useful constraints on cosmological parameters. In this 
paper we bring together the largest quasar samples avail¬ 
able so far with both optical/UV and X-ray observations, 
in order to build a sample with broad luminosity and 
redshift coverage, and a large enough size to obtain a 
meaningful estimate of the cosmological parameters. 

This paper is organized as follows. In Section we 
will describe the method adopted to obtain the Hubble 
Diagram for quasars, and to estimate cosmological pa¬ 
rameters (details are given in Appendix]^. In Section 
we discuss the general criteria for sample selection (witn 
a complete discussion presented in Appendix [A|). In Sec¬ 
tion 1 ^ we analyze the Lx — Lyv in our sample, and we 
validate its use as a cosmological probe. In Section 
we present the main results: the Hubble Diagram for 
quasars, and the measurements of cosmological param¬ 
eter it provides. We also compare and complement our 
results with those from supernovae, in particular by com¬ 
paring the Hubble diagram of supernovae and quasars 
below z^l.4, and calculating the cosmological parame¬ 
ters from joint fits. Our method and results are discussed 
in Section while in Section we will explore possible 
future extensions of our work, based on the inclusion 
of more already available sources to our sample, on new 
dedicated observations, and on forthcoming new surveys. 
Conclusions are outlined in Section |8l 

For luminosity estimates we adopted a concordance A- 
cosmolo gy with i?n = 70 kms~ ^ Mpc“^, Dm = 0.3, Da = 
1 — Dm ( Komatsu et al.|[2009 ). 

2. METHOD 


Our method is based on the non-linear relation be¬ 
tween Lx and Lyy: 

log(Lx) = 7log(Luv)- (1) 


From the we above equation we obtain 

log(Lx) = $(Fj/v,Da) = /3'-F7log(Luv)+2(7-l)log(L)i), 

( 2 ) 

where /?' depends on the slope and intercept (i.e. /?' = 

/3 -|- (7 — 1) log( 47 r)), Fx and L/jv are measured at fixed 
rest-frame wavelengths, and is the luminosity dis¬ 
tance, which in a standard ACDM model (i.e. in a cos¬ 
mological model with a fixed cosmological constant A) is 
given by 


L>l(z, Dm,Da) = sinh a/d// 

V Die 


dz 


(3) 


Jo FIo'\/Dm(I + Da -I- Die(l -I- 


where Die = 1 ~ Dm — Da. 

Our analysis will consist of two main parts: 

1) Validation of the relation at different redshifts. 
It is clear that any use of the above equations to deter¬ 
mine the cosmological parameters is based on the as¬ 
sumption of no evolution of the Lx-Lyy relation. While 
it is not possible to check whether a redshift dependence 
of the scaling parameter /3 is present (unless we have a 
physical model to determine it independently of the ob¬ 
served datm which is currently not the case, as discussed 
in Section^, we can test the linear shape of the cor¬ 
relation ana its slope 7 at different redshifts. In order 
to do so, we note that if Equation is analyzed in a 
sufficiently narrow redshift interval, the term containing 
Dl will provide a nearly constant contribution. In par¬ 
ticular, if its range of values within the chosen redshift 
interval is smaller than the intrinsic dispersion of the 
correlation, we will be able to test the relation by replac¬ 
ing luminosities with observed fluxes. An analysis of the 
function Dl{z) for several values of Dm and Da shows 
that in every redshift interval at z > 0.3 with amplitude 
A[logz] we have ApogDy] < 0.7A[logz]. 

The measured dispersion of the global Lx-Lyy in the 
main published works cited in the Introduction is of the 
order of 6 ~0.35-0.40. In order to have a negligible con¬ 
tribution by the distance term in a flux-flux relation, 
we need ApogDy] < 0.10 — 0.15. This implies a maxi¬ 
mum size of the redshift bins A[logz] < 0.1. This sets a 
first important requirement for a quasar sample in order 
to be well suited for cosmological studies: it should con¬ 
tain enough sources to allow a significant test of the Lx- 
Lyy relation in redshift bins smaller than A[logz] ~0.1. 

2) Determination of the cosmological parameters. 
We can derive the cosmological parameters from the re¬ 
lation in Equation!^ in two equivalent ways: 

• We fit Equation!^ by minimizing a likelihood func¬ 
tion (LF) consisting of a modified function, al¬ 
lowing for an intrinsic dispersion S: 


In(LF) 


^ I [log(J^x)i - ^{Fuv, DL)i]'^ 


(4) 

where $(Fyy,L>L)i is given by Eq. 2, and sf = 
tj/ -|- 5^, with (Ti and 6 indicating the measurement 
errors over Lx and the global intrinsic dispersion, 
respectively. We note that the dispersion 6 is 
much higher than typical values of try hence the 
fit is almost insensitive to the exact ai value. The 
results do not change if an additional error due to 
systematic effects is added to at, provided it does 
not exceed the intrinsic dispersion. In this case the 
free parameters are 6, the slope and intercept 7 and 
/3, and the cosmological parameters Dyf and Da 
( we assume no evolution of the equation of state 
of dark energy, i.e. w = —1). We note that the 
Hubble constant Hq is absorbed into the parameter 
/3: without an independent determination of this 
parameter, our fits are insensitive to the value of 
Ho- 


• We can use the best fit value and the uncertainty of 
the slope 7 (estimated in narrow redshift bins, as 
described above) to directly compute the distance 
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modulus for each quasar in the sample, simply re¬ 
arranging Equation in order to obtain the lumi¬ 
nosity distance as a function of the observed X-ray 
and UV fluxes. We then Ht Dl{^m,^a) as a. func¬ 
tion of redshift, minimizing a likelihood analogous 
to the one in Equation]^ In this case the free pa¬ 
rameters are the intrinsic dispersion, the cosmolog¬ 
ical parameters and Ha, and a residual scaling 
parameter /3'. 


The two methods outlined above are equivalent and pro¬ 
vide fully consistent results. However, we note that a di¬ 
rect fit of the luminosity distance as a function of redshift 
has two practical advantages: its easier visualization, be¬ 
ing formally a fit of a function of one variable, while the 
previous method requires a fit of Ex as a function of 
{F\jy,z), and its homogeneity with other cosmological 
probes using Hubble Diagrams (chiefly those based on 
supernovae la). 

The minimization has been performed using a Monte 
Carlo Markov Chain (MCMC) , with the emcee package 
(Eoreman-Mackey et al. 20131 in Python 2.7. Elat pri- 
ors were used for each parameter. In order to further 
check the results, we also fitted the data in two other 
ways: (1) we used a standard likelihood maximization, 
based on the Levenberg-Marquardt (LM) algorithm; (2) 
we performed a “brute force” maximization calculating 
the likelihood function in a grid of Hm, Ha, and <5 (the 
7 and j5 values minimizing LF can be computed analyt¬ 
ically as a function of the other three parameters). In all 
cases, we find fully consistent results. In the latter ap¬ 
proach, the confidence contours have been estimated via 
a bootstrap procedure. While the grid method is not effi¬ 
cient for precise parameter estimates (and has been used 
only to check the shape of the likelihood function), both 
the MCMC and LM methods have advantages (in par¬ 
ticular, the latter is faster). Here we prefer to report the 
results from the MCMC approach for consistency with 
the main published works on observational cosmology. 


3. THE SAMPLE 

To obtain an adequate coverage of both the Fyv — -2 
and Ex — z planes we need to combine wide, narrow, and 
deep-field surveys. In this first work we concentrated on 
data already available in the literature from published 
works on the aox ~ Euv correlation, with the require¬ 
ment of a good photometric coverage in the optical-UV, 
in order to obtain precise estimates of Fjjv and of the 
possible extinction due to dust. We sta r ted w ith the 
quasar sample presented by Steffen et al. (2006), which 
is a collection of several surveys such as the wide-field 
Sloan Digital Sky Survey (SDSS) quasar sample, low- 
redshift {z < 0.2) Seyfert 1 galaxies, and high-redshift 
{z > 4) optically selected AGNs. We ha ve expanded this 
sampl e with the quasars at z > 4 by [Shemmer et al.| 
(2006), the optically s e lected sam ple of high luminou s 
quasars by Just et al. (2007) and Young et al. (2010). 
These optically selected samples are then combined with 
the X-ray selected one analysed by iLusso et al .| (|2010[ 
for a total of ^1,138 sources (see appendixes I A. 1 1 aim 
for further details). 

All these catalogs provide a measurement of the rest- 
frame X-ray flux at 2 keV (Ex), and the UV flux at 
2500A (Euv)- In our work, we adopted the published 


AT 


X-ray fluxes for most sources, while we updated the Euv 
estimates by using whenever possible the available multi¬ 
color information compiled in the COSMOS, SDSS-DR7, 
and BOSS-DRIO quasar catalogs. While we refer the 
reader to Appendix for details on our sample, flux 
estimates and possible systematics, a brief summary is 
given below. 

The adopted catalogs include multi-wavelength data 
from mid-infrared to ultraviolet: MIPS 24 /im G03 
data, IRAC flux densities, near-infrared Y-J-H-K-bands 
(2MASS and/or UKIDSS), optical photometry (e.g. 
SDSS, Subaru, CFHT), and near- and far-UV bands 
(GALEX). Observed magnitudes are converted into 
fluxes and corrected for Galactic reddening by employ¬ 
ing a selective attenuation of the stellar continuum with 
Ry = 3.1. Galactic extinction is estimated individually 
for each object in all catalogs. Eor each source we con¬ 
sidered the flux and corresponding effective frequency in 
each of the available bands. The data for the SED com¬ 
putation from mid-infrared to UV (upper limits are not 
considered) were then blueshifted to the rest-frame and 
no K-correction has been applied. We determine a “first 
order” SED by using a first order polynomial function 
(i.e. straight line in the logv — log(izFj,) plane), which 
allows us to build densely sampled SEDs at all frequen¬ 
cies. This choice is motivated by the fact that a single in¬ 
terpolation with a high-order polynomial function could 
introduce spurious features in the final SED. Ejjv are ex¬ 
tracted from the rest-frame SEDs in the log iz — log(;/F;^) 
plane. In the case data do not cover 2500 A, fluxes are 
extrapolated at lower (higher) frequencies by consider¬ 
ing, at least, the last (first) two photometric data points. 

Finally, we corrected the E/jv estimates taking into 
account the redshift-dependent contribution of emission 
lines t o the photometric points, as discussed in Ap¬ 
pendix |A.4| An analogous correction to the values of 
Ex based on the comparison between our values and the 
results of a c ompl ete X-ray spectral analysis is discussed 
in Appendix |A.5| 


3.1. Source selection 

The whole group of quasars collected as described 
above is not yet well suited for cosmological studies. 
The most obvious sources deviating from the Ex-Euv 
relation are the heavily obscured Broad Absorption Line 
(BAL) and radio-loud quasars which have an additional, 
jet-linked. X-ray component. These sources were already 
removed by most of the parent samples, with a few ex¬ 
ceptions (mostly radio-loud objects) which we further ex¬ 
cluded from our sample. A few undetected BAL quasars 
may still be present at low redshift. These oWects may 
introduce a small bias in the Ex-Euv relatioi|/J which is 
at present unavoidable, unless UV spectra are obtained 
for the whole sample. 

The most serious issue potentially affecting a large frac¬ 
tion of the quasar sample is dust extinction. This effect 
is expected to be stronger at 2500 A than at 2 keV. As¬ 
suming a Galactic dust-to-gas ratio, a relatively small 
amount of dust, associated with a column density of a 

^ BALs are known to be X-ray obscured (e.g.,|Green et al.|l995| 
[Gallagher et all]|1999[ [Brandt et al.|[2000^ , and are not includea 
in previous studies ot optically selected samples because they can 
cause an artificial steepening of the Lx-L\jy correlation. 
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Fig. 1.— Distribution of the whole quasar sample (excluding 
radio-loud and BAL sources) in a Fi — F 2 plot, where Fi and F 2 
are the slopes of a power law in the \ogu — log(uFj^) plane, in the 
0.3-0.5 /i. and 0.3-1 fim intervals, respectively. The typical error 
for Fi and F 2 is of the order of 0.1. The red star represen ts the 
intrinsic quasar SED as estimated by (Richards et al.| (|2006|l. The 
solid red line is obtained by a ssuming increasing dusx extinction 
following the extinction law of|Prevot et al.|||1984|| and therefore, 
even if it nicely reproduces the correlation among the two colors, 
is not a fit to the points. The dashed lines, orthogonal to the solid 
one, show the colors corresponding to different values of the ratio 
rijY between the intrinsic and observed flux at 2500 A. 


few 10^° cm is enough to decrease the observed flux 
at 2500 A by a factor of ^two, while the flux at 2 keV is 
nearly unaffected. 

Since dust extinction is expected to redden the ob¬ 
served optical-UV spectrum, in order to analyze this 
issue we computed for each object the slope Ti of a 
\og{i') power law in the 0.3-1 /im (rest frame) 

range, and the analogous slope r 2 in the 0.3-0.5 /im 
range (rest frame). The Ti — r 2 distribution is shown 
in Figured As expected, the two slopes are in general 
well correlated. The large tail of the Ti — r 2 distribution 
towards low values (i.e. redder spectra) is indicative of 


iPrevot et al.||1984 

appropriate for un- 

[Hopkins et al. 

20U4 

Salvato et al.||2009l 

i 'l — I 2 correi 

ation as a function of ex- 


tinction. We obtained the red solid line shown in Fig¬ 
ure where the starting point (with zero e xtinction, 


plotted with a star) is derived from the SED of Richards 
et al.| ( |2006[ i.e. Fi = 0.95, F 2 = 1.35). This line is 


consistent with the best fit line to the observed points, 
and confirms that the main physical driver of the large 
color distribution is dust reddening. This implies that 


the observed fluxes at 2500 A are in many cases signifi¬ 
cantly underestimated. The dashed lines, orthogonal to 
the solid one, show the colors corresponding to different 
values of the ratio ruv between the intrinsic and ob¬ 
served flux at 2500 A. The corresponding attenuation at 
2 keV, assuming a Galactic dust-to-gas ratio, is instead 
negligible (about 5% if the ratio between intrinsic and 



Fig. 2.— Redshift distribution of the quasar sample used for the 
cosmological analysis. The total number of sources is 808. 



log Tuv [erg s ^ cm ^ Hz 


Fig. 3.— Fuv ~ -^x distribution for our sample, with colors and 
point shapes representing different redshift intervals, as indicated 
in the plot. The correlation between .F’uv and Fx is due to the 
physical correlation between the corresponding luminosities. The 
observed, large dispersion is due both to the intrinsic dispersion of 
the I/uv ~ relation and to the redshift distribution at each flux 
interval. The information on the cosmological parameters is there¬ 
fore encoded in this dispersion, as visually shown by the three red¬ 
shift intervals, with the lowest (blue diamonds) and highest (green 
squares) redshift objects lying preferentially on the upper and lower 
envelopes of the distribution, respectively. 


observed flux at 2500 A is ruv = 10). 

In order to obtain an optimal sample for the analysis 
of Fx-F\jv relation, we need to correct the observed UV 
fluxes for dust extinction. This correction likely intro¬ 
duces further dispersion to the correlation, due to the 
spread in the intrinsic SED of quasars, and the possi¬ 
ble differences in the dust properties (and hence in the 
extinction law) among quasars. Such uncertainties are 
likely to be larger when the required extinction correc- 
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tion is also larger. We therefore need to choose an opti¬ 
mal balance between including more objects (so increas¬ 
ing the statistics) and applying smaller extinction correc¬ 
tions (so reducing the dispersion). This choice does not 
introduce biases in the subsequent analysis, and can be 
made following the criterion of the smallest errors in the 
final best fit parameters of our analysis. We therefore re¬ 
peated the cosmological analysis (discussed in detail the 
next Section) for several choices of the maximum allowed 
rl;v^ and we selected the sample with the smallest errors 
in the determination of and Ha (specifically, the one 
minimizing the area of the l-tr Hjvf — Ha contour). As 
noted above, absorption effects on the X-ray fluxes are 
expected to be less important. However, some extreme 
case are possible, for example if a very low dust to gas 
ratio implies low UV extinction and relatively high X- 
ray gas absorption. In these cases the spectrum is much 
flatter than in standard quasars, with photon indexes 
r < 1 . At the other extreme, in a few cases it is possible 
to observe very steep spectra T > 3 due to contami¬ 
nation by galactic sources in galaxies with a very high 
star formation. In order to exclude these extreme out¬ 
liers, we required the photon index of our sources to be 
in the standard range for unobscured, luminous quasars, 

r = 1.8± 1. 

The final “best sample” is summarised in Table and 
it consists of 808 objects with ruv < 10 (BAL and radio- 
loud quasars have been neglected). In Figures and 
we plot the redshift and Fuv ~ -Px distributions. Opti¬ 
cal and X-ray luminosities/fluxes with their uncertain¬ 
ties are listed in Tablej^ (different subsamples are flagged 
from 1 to 5, see Note m Table [^. 

4. ANALYSIS OF THE Lx - Lyy RELATION 

Based on the method outlined in Section we divided 
our sample in narrow redshift bins, in order to check the 
redshift dependence of the Lyv ~ Ax relation. The two 
requirements for such analysis are ( 1 ) that the scatter due 
to the different luminosity distances within each bin is 
small compared with the intrinsic dispersion and ( 2 ) that 
each bin is sufficiently populated for a meaningful check. 
The first condition is met for bins equally spaced in log z, 
with Alogz <0.1. The second condition is fulfilled for 
redshifts z > 0.5, while for lower redshifts (and narrower 
intervals, due to the first condition) we would have less 
than 15 objects in each bin (Figure 2). Therefore, we 
performed this analysis only in the z = 0.5 — 6.5 range. 

We divided the sample in 12 intervals with Alogz = 
0 . 1 , and for each redshift interval we performed a linear 
fit of log Fx - log Fuv relation, log Fx = 7 ^ log Fyy + Pz, 
with free and Pz- The results for the parameter jz are 
shown in Figured where we plot the 7 z — z relation, and 
three examples of our fits in three intervals centered at 
redshift z ~1,2 and 4. The value of jz is consistent with 
a constant value at all redshifts, without any significant 
evolution. The average gives 7 = ( 7 ^) = 0.60 ± 0.02. 

5. DETERMINATION OF THE COSMOLOGICAL 
PARAMETERS 

We adopted the average value of 7 obtained in the 
previous Section, 7=0.60±0.02, to estimate a distance 
modulus, DM = 51og[(FA/(10 pc)] for each quasar in 






Z 


Fig. 4.— Fx-Fuv correlation in narrow redshift intervals. The 
dispersion due to the different distances within each redshift in¬ 
terval is negligible with respect to the intrinsic dispersion of the 
Fx — F\jy correlation (see text for details). Therefore, the relation 
between fluxes is a good proxy of that between luminosities. Up¬ 
per panel: examples of the correlation in three redshift intervals 
at 2 ~1,2 and 4, respectively. Lower panel: best fit values and 
dispersion S of the correlation slope, 72 , in each redshift interval. 
The horizontal lines show the average values, i.e. ( 72 ) = 0.60di0.02 
and {6) = 0.3. 

our sample, defined from Equations 2, 3 as: 

DM = ^^ [log(Fx) - 7log(Fuv) - P']- (5) 

Since the intrinsic value of P in not known, P' can be 
considered an arbitrary scaling factor. 

In Figure we plot our results for each quasar, and 
for averages in narrow redshift bins (for visual purpose 
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TABLE 1 

Summary of the best quasar sample. 


Number 


Z 

log Luv 
erg sec“^ Hz“^ 

log Lx 

erg sec“^ Hz“^ 

Reference 


222 

0.061 

- 6.280 

28.94 - 32.90 

25.14 - 28.09 

Steffen et al. 

2006 

1 

20 

1.760 

- 4.610 

32.14 - 32.85 

27.14 - 28.07 

Just et al. ( 



14 

4.720 

- 6.220 

31.38 - 32.42 

26.62 - 27.62 

S 

lemmer et a 

.Nii 0 o 6 


327 

0.345 

- 4.255 

28.96 - 31.79 

25.08 - 27.14 

Lusso et al. 

2010 


225 

0.173 

- 4.441 

29.33 - 32.40 

25.08 - 27.78 

Young et al. 

o 

o 

1 

808 

0.061 

- 6.280 

28.94 - 32.90 

25.08 - 28.09 

Total 




TABLE 2 

Optical and X-ray properties of the “best sample” . 


Name 

RA 

DEC 

Redshift 

log Luv 
[erg s~^Hz“^] 

log Lx 

[erg s~^Hz“^] 

log Luv 

[ergs“^cm~^Hz“^] 

log Lx 

[erg s~^cm~^Hz“^] 

Group'' 

022356.30-085707.8 

35.985 

-8.952 

1.575 

30.92 

26.67 

-27.29 ± 0.01 

-31.54 ± 0.10 

1 

022435.92-090001.4 

36.150 

-9.000 

1.611 

30.97 

26.76 

-27.27 ± 0.01 

-31.47 ± 0.11 

1 

023306.25-1004614.5 

38.276 

0.771 

2.292 

30.66 

26.96 

-27.94 ± 0.03 

-31.65 ± 0.07 

1 


Notes—This table is presented in its entirety in the electronic edition; a portion is shown here for guidance. 

^ Group flag: [Steffen et al.| ( [2006| l (l), |Just et al.| ( |2007[ l (2), jShemmer et al.| ( |200^ (3), XMM—COSMOS (4), [Young et al^ | |2010| | (5). 



z 

Fig. 5.— Hubble Diagram for the quasar sample (small gray points) and supernovae (cyan points) from the Union 2.1 sample (Suzuki et 
al. 2012). The large red points are quasar averages in small redshift bins. The inner box shows a zoom of the z = 0 — 1.5 range, in order 
to better visualize the match between the SNe and the quasar samples. The continuous line is obtained from a joint fit of the two samples 
assuming a standard ACDM cosmological model. The relative normalization between SNe and quasars is a free parameter of the fit, and 
it is estimated with an uncertainty lower than 1%. 
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For our joint quasar-SNe joint analysis we used the 
Union 2.1 SNe sample of the Supernovae Cosmology 
projecl]^ (Suzuki et al. 2012). We derived the cosmo¬ 
logical parameters from the Hubble diagram of SNe only 
(obtaining the same results as in Suzuki et al. 2012, and 
so demonstrating the full consistency of our analysis with 
the published one) and for the joint Hubble Diagram. 
The results for the determination of the cosmological pa¬ 
rameters and Da are shown in Figure The 68 % 
and 95% contours for Dm and Da, assuming a standard 
ACDM model, are derived from our analysis of the Hub¬ 
ble diagram of quasar s only (blue contou rs), from the 
SNe Union 2.1 sample Suzuki et al. (2012 empty black 
contours), and from a joint ht quasars -|- SiNe (orange-red 
contours). Thanks to the large redshift overlap between 
the SNe and quasars samples (inner box of Figure 5) the 
cross-calibration parameter /3’ is estimated with an un¬ 
certainty lower than 1%. We note that the contour plot 
relative to the joint sample is not a statistical intersec¬ 
tion of the two individual contours, but the result of a 
simultaneous fit of the DM — z relation for the combined 
sample. 


6 . DISCUSSION 


Fig. 6.— 68% and 95% contours for Qm and £7 ai assuming a 
standard ACDM model, as derived from our analysis of the Hubble 
diagram of q uasars (blue), from the SNe Union 2.1 sample |Suzuki| 
|et al.| | 2012 [ | (empty black), and from a joint fit (orange-red^ 

only). We fitted the DM — z relation as discussed in 
Section]^ assuming a ACDM model, with Da, Dm, tbe 
scaling parameter /?' and an intrinsic dispersion S as free 
parameters. We also performed a direct fit of Fx as a 
function of Fuv and z, with 7 , (3, Da, Dm, and 5 as free 
parameters. The two methods provide fully consistent 
results (as expected, since they are two different rear¬ 
rangements of the same data). A summary of the results 
of cosmological fits is provided in Table [3 The best fit 
value of the dispersion is (5=0.30, lower man the values 
previously reported in the literature. This is a positive 
effect of our source selection, which removed from the 
sample many outliers with poor quality data and/or af¬ 
fected by absorption. 

One of the most interesting aspects of the analysis 
through the distance DM — z relation is that it provides 
a direct way to compare and merge our data with those 
from supernovae (SNe). We can of course compare our 
constraints with those from any cosmological measure¬ 
ment, such as CMB, BAO, clusters, etc. However, this 
analysis can be done only on the final values of the cosmo¬ 
logical parameters obtained from each method indepen¬ 
dently, since each method is based on the measurement 
of different physical quantities (CMB fluctuations, bary- 
onic oscillations, etc.). In our case instead, both SNe and 
quasars are used as “standardized candles” and the cos¬ 
mological parameters are obtained from the same physi¬ 
cal quantity, i.e. the distance modulus. Technically, the 
only difference is the absolute calibration, or the “zero- 
point” of the DM — z relation, which can be measured 
for SNe, but is unknown for quasars. As a consequence, 
a joint ht of the DM — z relation for SNe and quasars 
requires that the j3' parameter in Equation m is used as 
a free parameter of the ht in order to cross-cSlibrate the 
DM of the two group of sources. 


The main results of our work are the validation of the 
linear log Lyv — log Lx relation at all redshifts, and its 
application as a cosmological probe. The physical origin 
of the relation is unknown. Here we only notice three 
relevant aspects. 


1. The analysis in narrow redshift bins has shown no 
signihcant redshift evolution of the shape of the 
relation (a linear ht is adequate at all redshifts) 
and of the correlation slope, 7 . We cannot di¬ 
rectly test the constancy of the parameter /5, but 
an important consistency check of this point is pro¬ 
vided by the comparison between the Hubble dia¬ 
grams of quasars and SNe at z < 1.4 (the inner 
panel in Figure [5). The perfect match between the 
two curves is indicative that the evolution of the 
/3 parameter is negligible. Either that, or there 
is some unknown redshift evolution is also present 
between luminosity and redshift in SNe, with the 
same dependence as for the logLuv — logTx for 
quasars. The need to invoke this “physical con¬ 
spiracy” is enough to rule out this possibility. For¬ 
mally, the possibility of a redshift evolution of the 
logLuv ” log Ax relation beyond z = 1.4 remains, 
even if it appears quite unlikely, and not based on 
any physical motivation. On this regard, we notice 
that despite a large dispersion in many observa¬ 
tional properties of quasars, their optical-UV SED 
does not show any evolution with redshift and/or 
luminosity, up to the most extreme cases, as shown, 
for example by the similarity between the average 


quasar SED ([Richards et al.|2006 Hao et al. 


and the SEP of a luminous z ^ 7 quasar ([iVlort- 


lock et al. 2012). The same absence of evolution¬ 


ary/luminosity effects is observed in the s pectral 
properties of quasars in har(l X-rays (e.g. Young 


et al. | |2010[ [Lusso et al.||2010[ ). 


^ http://supernova.lbl.gov/Union/ 
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TABLE 3 

Results of cosmological fits. 



ACDM® 

juCDM^ 


tDQ'UIaCDM'^ 



SIa 

fljvf 

W 


WQ Wa 

DATA*" 


^•^^-0.30 





DATA'^+SNe 


n 'To + O.OS 





10 ® QSO'^ 

fj q'7+0.02 
^•^'-0.02 

^•oy_0.07 





10® QSO'^-l-SNe 

n 97’+0.02 
'^•^'-0.02 

'^-0.04 

0 26'*'^'^^ 
'-'•^0-0.014 

'4-yO_0.05 



10 ® QSO*" 

n o'7+0.02 
^•^'-0.01 

0 73+0-01 

A Q7+0 02 

-1 

A Q’T+U.Ud 
-0.02 

1 AO+U.Ub q + U.4 

-L-^o 0.05 ^•'^-0.4 


Note. — a: Data samples: DATA: our actual sample; 10® QSO: the simulated 10® QSO sample from the cross correlation between the 
SDSS DR7 and DRIO quasar catalogs and the eROSITA all-sky survey; 10® QSO: simulated sample of one million quasars from future all 
sky surveys (see text for details). 6 : Standard ACDM model, c: Cosmological model assuming a flat Universe and a non-evolving dark 
energy parameter w. d: Cosmological model assuming a flat Universe and an evolving dark energy equation of state: w = uiq + Wa x (l-l-z). 


2. The cosmological application presented in this pa¬ 
per is not directly based on the physical interpre¬ 
tation of the Luv ~ Lx relation, but only on its ob¬ 
servational evidence. This situation is qualitatively 
similar to what is presently known about super¬ 
novae as cosmological tools: the “standardization” 
of the supernovae luminosity is obtained through 
the empirical relation between the peak luminosity 
and t he slope of the luminosity light curve after the 
peak ( Phillips! 1993 1. 


3. The large observed dispersion in the Luv — Lx re¬ 
lation (~ 0.30 in logarithmic units, i.e. a factor of 
two in physical units) is obviously the main limi¬ 
tation for precise cosmological measurements. We 
note for example that the dispersion in the distance 
modulus-redshift relation is much smaller for su¬ 
pernovae than for quasars. In other words, quasars 
as standard candles are much less precise than su¬ 
pernovae. However, the main advantages in using 
quasars are the wider redshift range, and the pos¬ 
sibility of future large improvements of the sample. 
This is discussed in detail in the next Section. 

The origin of the intrinsic dispersion must be related to 
the physical origin of the relation, and is therefore un¬ 
known. However, the observed dispersion may be sig¬ 
nificantly larger than the intrinsic one. Here we briefly 
discuss possible contributions to the observed dispersion: 

- Inclination effects: assuming a disk-like UV emission 
and an isotropic X-ray emission, we expect that edge-on 
objects have a relatively low o bserved UV-X-ray lumi¬ 
nosity ratio (Risaliti et al. poTTj ). This implies that edge- 
on quasars are outliers in the logLuv ~ logTx relation, 
and a contribution to the dispersion comes from the dis¬ 
tribution of disk inclinations. In an optical-UV selected 
sample the distribution of inclination angles 9 has a peak 
value at cos6> = 5/9 , with a dispersion a{cos9) = 2/9 
(Risaliti et al. 2011). On average, this ^40% contribu- 
tion to the dispersion is negligible with respect to the to¬ 
tal one (after correction, the observed dispersion would 
decrease from 0.30 to ^0.29), but this effect may be re¬ 
sponsible for some of the largest outliers. In an X-ray 
selected sample, there would be no selection bias with 
respect to the disk inclination (assuming isotropic X-ray 
emission), and the dispersion due to the inclination effect 
would be slightly larger, but still negligible with respect 
to the total one. 


- Variability: our optical-UV and X-ray observations 
are not simultaneous, therefore variability must add a 
significant contr ibution to the o bserv ed dispersion. Pre¬ 
vious studies of Gibson et al. (2008) on SDSS quasars 
suggest that variability cannot explain all of the ob¬ 
served dispersion in the Oox — Duv relation. It is how¬ 
ever possible that simultaneous observations (for exam¬ 
ple, taking advantage of the Optical Monitor on-board 
XMM-Newton, and/or the XRT and UVOT instruments 
on-board Sw ift) could reduce the intrins ic dispersion, as 
suggested by Vagnetti et al. (2010 2013). We note how¬ 
ever that even if variability is responsible for a significant 
fraction of the observed dispersion, simultaneous obser¬ 
vations may be not enough to obtain a tighter correla¬ 
tion, if the delays between optical and X-ray variations 
are longer than the typical observing time. This seems 
to be the case in low -luminosity AGN, such as NGG 5548 
(Edelson et al.|2015 ) where the X-ray emission leads the 
optical one by about a day. It is reasonable that in higher 
luminosity quasars the delays are significantly longer. 

- Absorption correction. We corrected our sample as 
explained in Section 3, assuming an intrinsic quasar SED 
as in Richards et al. (2006), and a standard extinction 
law. This simple approach has several limitations: 

1) We neglected possible contributions from the host 
galaxy. While this is likely a safe assumption for the 
X-ray emission, it is possible that the contribution of 
the galaxy in the optical-UV cannot be ignored, espe¬ 
cially for low luminous AGN {Lbol < 10"^^ erg s“^). 
This would alter both the estimate of the UV flux of 
the quasar, and the absorption correction. We ex¬ 
pect this contribution not to be relevant for most ob¬ 
jects, given the high luminosity of our sources (aver¬ 
age log(LBOL)=46. a, and the good match between the 
Ti — r 2 distribution with the expectations from dust ex¬ 
tinction (Eigure 1). 

(2) We assumed the same SED for all sources, while 
intrinsic differences are present even in blue color/UV- 
selected samples (Elvis et al. 1994, Richards et al. 2006). 

(3) We assumed the same dust extinction law for all 
sources, neglecting possible differences in the composi¬ 
tion of dust grains in each object. 

All these simplifications likely contribute to the observed 
dispersion, which may be reduced with a future, more 
detailed analysis of the optical-UV SED of each source. 

® Estimated from the luminosity at 2500 A and assuming a bolo- 
metric correction of 3, as in Krawczyk et al. (2013). 




















A Hubble diagram for quasars 


9 


- Systematic effects. A final concern regarding our find¬ 
ings is related to possible systematic effects in the sam¬ 
ple selection. Specifically, the fits of the Hubble diagram 
may be biased in different ways if the observations in one 
of the two bands are biased towards high values. In our 
sample, such a bias may be present due to our choice to 
include only sources with both UV and X-ray measure¬ 
ments, avoiding upper limits. The Steffen et al. (2006) 
and Young et al. (2010) subsamples, both optically se¬ 
lected, contain a small fraction of X-ray non-detections 
(12% and 13%, respectively). Analogously, the Lusso et 
al. (2010) X-ray selected subsample could be undersam¬ 
pled at low UV fluxes. In the latter case the effect on the 
fit of the Fx-F[/v relation is expected to be smaller than 
for UV-selected samples, due to the broader flux range 
in the optical-UV than in X-rays (Figure [2). In order to 
further investigate this issue, we simulate the effect of 
cutting the available sample at a relatively high X-ray 
flux values (log(F’x,min) = —32 in the same units as in 
Figure]^. This is equivalent to assuming that our sample 
is derived from the cross correlation between an optically 
selected catalog and an X-ray survey with a uniform X- 
ray lower limit Fx^min- The results are a significantly 
lower slope of the logLyv — log Ax relation (7 ^ 0.5), 
and a distorted Hubble diagram, providing cosmological 
parameters biased towards lower values of both D,m and 
Ha. Making a general statement from the result of this 
test is not straightforward, since the effects of such biases 
on the estimates of the cosmological parameters depend 
on both the redshift and flux distributions of the sam¬ 
ple. Different redshift distributions may have a different 
impact of the flux limit on the Hubble diagram. It is 
also difficult to quantitatively estimate the possible ef¬ 
fects of these kinds of biases on our findings, given the 
non-homogeneous composition of our sample. However, 
we are confident about our results, based on the following 
two considerations: 

1. The results of the analysis of the '^ox-^uv 
relation in Young et al. (2010) are insensitive to 
the inclusion/exclusion of upper limits, suggesting 
that the bias is small. 


2. We checked for the possible effects of these bi¬ 
ases by cutting our sample at different UV and X- 
ray minimum fluxes, starting from log(Fx,min) = 
—33.5 and log(Fuv,min) = —29. We already 
discussed the case with a cut at log(Fx,min) = 
—32, where the effect of the bias was not negli¬ 
gible. However, we verified that the results are 
not significantly altered up to minimum fluxes of 
log(Fx,min) ~ —32.5. This shows that our current 
sample is not seriously affected by systematics re¬ 
lated to flux limits. We note that this may become 
a much more serious issue with future, larger sam¬ 
ples which will allow more precise measurements of 
the cosmological parameters, and will therefore be 
sensitive to smaller systematic effects. 

7. FUTURE DEVELOPMENTS 

The analysis presented in this paper already provided 
important results: the validation of the Auv — Ax 
relation at all redshifts, and the extension of the Hubble 


Diagram up to z ~ 6 . However, the constraints on the 
cosmological parameters are still loose compared with 
other methods, and the improvement on the measure¬ 
ment errors obtained by combining our constraints to 
those from SNe, CMB, IIAOs and clusters is not signifi¬ 
cant. However, the present work is mainly intended as a 
demonstration of the power of the method, while most 
of its potential is still to be exploited in future work, 
both with available and forthcoming data. Here we list 
the main foreseeable future developments. 


7.1. Sample improvements - present. 

The work presented here is based on the data available 
in the literature, from samples already used to study the 
Q!ox — Auv correlation. We plan to extend the analysis 
in at least three directions. 

1. We will increase the sample size at low redshift. 
The present shape of the Hm — Ha contour, elon¬ 
gated along the Ha direction, is due to the poor 
sampling at low redshift (10% of the sample is at 
z < 0.5, only 6 objects are at z < 0.1), where the 
effect of Ha is higher. Our simulations show that 
an optimal redshift distribution of the quasar sam¬ 
ple should be about constant in logz. An increase 
of the number of sources at low redshift can be ob¬ 
tained by adding local bright Seyfert galaxies and 
quasars with X-ray and optical-UV observations. 
The crucial point for nearby and/or low luminosity 
sources, which requires a source by source check, 
is the accurate determination of the flux at 2500 A 
taking into account for possible galactic contamina¬ 
tion, dust extinction, aperture effects. We expect 
to add about 100 sources from different catalogs in 
the literature at z < 0.1. We note that this addi¬ 
tion will be important in order to improve the preci¬ 
sion in the estimate of the cosmological parameters 
with quasars alone, while the constraints from the 
total Hubble diagram of SNe and quasars will re¬ 
main unchanged. Indeed, the constraints from SNe 
are at present much more stringent in the overlap¬ 
ping part of the diagram (i.e. at z < 1.4). There¬ 
fore, the quasars at z < 1.4 are useful (a) to test 
the method and (b) to cross-calibrate the SNe and 
quasar diagrams (we recall that this is needed be¬ 
cause the scaling constant (3 of the Ax-Ayv cor¬ 
relation in quasars in unknown). Since our sam¬ 
ple already provides an excellent cross-calibration 
at z < 1.4, only new sources above z = 1.4 will 
contribute to increase the measurement precision 
if SNe and quasars are used together. 

2. We will analyze through simulations the effect 
of adding to the sample bright, intrinsically blue 
quasars at high redshift (z > 4) from the SDSS 
surveys, and propose them for new X-ray observa¬ 
tions. This will greatly improve the quality of the 
Hubble diagram in this redshift interval, and will 
provide a precise test of the cosmological model in 
a redshift range where no other method is available. 

3. We will cross-correlate the SDSS quasar catalogs 
with the XMM-Newton and Chandra archives. We 
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Qm Qm 

Fig. 7.— Simulations for a sample of 100,000 quasars obtained cross-correlating the SDSS DR7-I-DR10 quasar catalogs with the future 
eROSITA X-ray all-sky survey (see te xt for details) . Left p anel: 68% and 95% contours for Qm ^ ACDM model (blue: quasars; 

empty black: SNe Union 2.1 sample |Suzuki et al |(|2012^; orange-red: joint fit). Right panel: same for a model with the dark energy 
equation of state parameter w, assumed to be constant at all redshifts {w = —1 is equivalent to having a cosmological constant). 





Fig. 8 .— Simulations for a sample of 1,000,000 quasars with the same Fvv-z distribution as our present sample (blu e contours), compared 
with the b est measurements available today, obtained from Planck, lensing, BAO and supernovae measurements (|Planck Collaboration! 
|et al.|2015l empty red contours). Left panel: 68% and 95% contours for Qm and Qa in a ACDM model; Middle panel: same tor a model ot 
a tlat Universe with the dark energy equation of state parameter w, assumed to be constant at all redshifts (m = — 1 is equivalent to having 
a cosmological constant). Right panel: same for a model of a flat Universe with evolving dark energy equation of state, w = wq + (f + z)wa ■ 
The mismatch between the contours are due to the difference between the si mulated values (ACDM model with Qm = 0.27 and Qa = 0.73) 
and the best estimates from the Planck+BAO+lensing+SNe measurements ||Planck Collaboration et al.|20f^ Qai ~ 0.31, S2a ^0.69). 


expect to obtain several hundred sources with 
optical-UV and X-ray coverage, which will greatly 
improve our statistics. Within this project, we plan 
to compare the results obtained using the SDSS 
phot ometry with those from the Optical Moni¬ 
tor (Page et al. 20121 on-board XMM-Vewtoti, 
which are simultaneous with the X-ray observa¬ 
tions. This test will be important to investigate 
the effect of variability on the observed dispersion 
of the Luv ~ Lyi correlation. 

7.2. Sample improvements - near future. 


The final SDSS-III quasar catalog (Alam et al. 
contains at present almost 300,000 objects, h'or al 


20151 

these 


sources we can obtain a 2500 A flux, and a redden¬ 
ing/extinction estimate, as done for the present sample. 


A major increase in the number of X-ray observations 
of quasars will be crucial for the future use of quasars as 
cosmological probes. The eROSITA (extended ROe ntgen 
Survey with an Imaging Telescope Array) satellite (Mer- 


loni et al.||2012p is expected to survey the whole X-ray 

10“^^ erg s“^cm“^ (^2500 quasars over 
,2 


sky down to 
,2 


18 deg^, i.e. ^ 140/deg^). The fraction of Sloan quasars 
detected at the eROSITA limiting flux is ^ 25%: essen¬ 
tially 1/4 of all SDSS quasars are expected to be detected 
by eROSITA (Menzel et al. 2015, MNRAS submitted). 
We applied the same color selection as in Figure H to 
the whole SDSS DR7H-DR10 samples, and simulated an 
X-ray flux measurement for each object. We assumed 
the observed linear log Lx — log Luv relation, with slope 
7 = 0.60 and dispersion 6 = 0.30, and a ACDM cosmol- 
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ogy with = 0.27 and Ha = 0.73. Finally, we conser¬ 
vatively assumed that about half of the SDSS quasar will 
have an eROSITA X-ray flux, obtaining a final sample of 
^100,000 quasars (over the whole area covered by SDSS 
and eROSITA, i.e. 8,000 deg^). We then fitted to this 
sample a ACDM model, and a model with a dark energy 
equation of state parameter w (assuming a flat Universe). 
The results are shown in Figure and demonstrate the 
potential of quasars as cosmological probes: the preci¬ 
sion that can be achieved with SNe and quasars (i.e., 
only using a Hubble diagram of “standard candles”) is 
similar to that obtained today by co mbining all the avail- 


able method s (see, e.g. Figure 6 of Suzuki et al. 2012 
Figure 16 of|Betoule et al.||2014|). 

The simulation discussed here assumes only statistical er¬ 
rors, thus neglecting the possible systematic effects due 
to contribution of emission lines to Fuvi to extinction 
correction, and to the flux limits in the sample selection. 
We already discussed these points for the present sam¬ 
ple, concluding that their effects are negligible. However, 
the situation may be different with a much larger sample 
and more precise estimates of the cosmological param¬ 
eters. In order to further investigate these issues, we 
performed the following checks: 

1) We simulated the X-ray fluxes for the sample of SDSS- 
DR7 quasars of Shen et al. (2010), starting from their 
estimates of Fuv ■ We then replaced the Uyv values with 
the ones obtained applying our method based on pho¬ 
tometric points, and performed the cosmological analy¬ 
sis. We obtained cosmological parameters significantly 
shifted with respect to the ones assumed in the simula¬ 
tion {V,M ^ 0.15±0.03, to be compared to the simulated 
value of 0.27, while Ha is marginally compatible with 
the simulated value of 0.73 ). If a correction like the one 
discussed in Appendix |A.4| is applied, the discrepancy is 
reduced, but is still present. We conclude that, in order 
to achieve the precision shown in Figurej^ a careful spec¬ 
tral analysis is needed. This problem may be reduced by 
a more detailed SED fitting based on the available SDSS 
photometric points complemented by NIR (when avail¬ 
able) and WISE data. We will further investigate this 
point in forthcoming papers. However, in the specific 
case of the SDSS-beROSITA quasar sample, we will be 
able to precisely estimate Fuv ^^nd Fx from an analysis 
of the UV and X-ray spectra, which will be available for 
all the sources. An analogous check on the possible bias 
introduced by the extinction correction is less straight¬ 
forward, because it is not clear whether the uncertainty 
in such correction is redshift-dependent. Here we only 
notice that the statistics will be enough to restrict the 
sample to the bluer objects, and that a more precise cor¬ 
rection can be obtained from a complete spectral analysis 
rather than from the photometric colors. 

2) A similar issue as discussed above may be present 
in the determination of the X-ray fluxes. A simulation 
analogous to the previous one has been done based on 
the da ta in Fig. 9 (right panel, see discussion in Ap¬ 
pendix A.5|. In this case however we do not find any 
significant shift in the best fit values of the cosmological 
parameters. Furthermore, we expect to be able to obtain 
the value of Fx from a complete spectral analysis for all 
the sources in the sample. 

3) The flux limits of our sample may introduce a bias 


in the observed Lx-Ljjv relation: objects with an ex¬ 
pected (based on the relation) flux near to the sample 
flux limit will be observed only in case of positive fluctu¬ 
ations. This issue is expected to be relevant mostly in the 
X-rays, due to the larger flux range in the UV than in 
the X-rays (which in turn is due to the slope of the corre¬ 
lation being < 1). There are two ways to deal with this 
issue: the first one is to include non-detections in the 
analysis. This will require a different statistical treat¬ 
ment, in order to properly account for censored data. A 
second way to obtain an (almost) unbiased sample is to 
include only objects that would be observed even in case 
of negative flux fluctuations. To do so, for each Fuv we 
considered the expected 2 cr lower limit of Fx, based on 
the distribution in Fig. 3. We then included the object in 
the sample only if this minimum X-ray flux is above the 
detection limit (regardless of the observed value of Fx)- 
We checked the effects of this cut varying the rejection 
threshold, and we found that, considering the expected 
X-ray flux limit of the eROSITA all-sky survey, we may 
have a bias in the estimate of the correlation slope 7 , 
which in turn alters the shape of the Hubble Diagram. 
We found that a filter like the one presented above, with 
a 2(7 threshold, is enough to remove the bias. The frac¬ 
tion of rejected objects is of the order of 20%. Consid¬ 
ering that the total number of SDSS quasars is close to 
300,000, even applying all the cuts discussed above, the 
number of quasars in the final, clean sample will easily 
remain higher than the value of 10 ® we assumed in the 
cosmological simulations. 


7.3. Sample improvements - in 10-15 years. 

We can extend our simulations trying to predict the 
available data provided by the next generation of ma¬ 
jor observatories. For example, the ESA mission Euclid 
will identify a fe w million quasars th rough slitless near- 
IR spectroscopy ( Laureijs et al.|20r2 ), with thousands at 
very high redshifts (z > 6 ) ; the Large Synoptic Survey 
Telescope (|Ivezic et al. 120081 LSST) will discover through 
photometry and variability millions of quasars, with an 
efficiency and c ompleteness analogous to that of spectro¬ 
scopic surveys (|Edelson &: Malkan||2012|). On the X-ray 
side, the Athena (Advanced 'Feiesciope tor High ENergy 
Astrophysics) observatory will have a large held of view 
imager, with resolution of a few arcsec, capable of detect - 
ing unabsorbed quasars up to z ~ 8 (Aird et al. 2013). 
It is expected that a significant fraction of the observing 
time of Athena will be used to perform wide area sur¬ 
veys. The match between these two observatories will 
likely provide samples of at least several hundred thou¬ 
sands quasars with both optical/UV and X-ray measure¬ 
ments. Even better. X-ray su rvey telescopes, such as the 
Wide-Field X-ray Telescope ( Murray et al.||2009 ) would 
further increase these numbers by an order cif magnitude: 
about 10^ quasars are expec ted in the WFXT surveys, 
with ^1,600 objectw at z >6 (Gilli et al.||2010). In order 
to have a first hint of the possible use of such samples as 
a cosmological estimator, we simulated a sample of a mil¬ 
lion quasars, with the same properties as the presently 
available one. The results are illustrated in Figure]^ and 
show how the Hubble diagram of quasars may become a 
fundamental tool for precision cosmology, and in partic¬ 
ular for the determination of possible deviations of the 
dark energy component from the standard cosmological 
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constant. We note that this simulation is optimistic in 
the sense of assuming the availability of these new ob¬ 
servational facilities, but is quite conservative in other 
respects. A million quasars in probably a large underes¬ 
timate if Euclid, LSST, and a large area X-ray surveyor 
all become available. The simulated redshift distribution 
(obtained from the presently available sample) does not 
include the thousands of quasars at redshift z > 6 ex¬ 
pected from Euclid, which would in particular improve 
the precision of the determination of m (which, in turn, 
is partly degenerate with the Wa parameter in models 
with an evolving equation of state of the dark energy, 
W = Wo + Wa X z/{1 + z)). 


8 . CONCLUSIONS 

We have demonstrated that the non-linear relation be¬ 
tween UV and X-ray luminosity in quasars can be used 
as a cosmological probe. We have shown that within 
the precision allowed by present data, this relation does 
not show any evolution with redshift and/or luminos¬ 
ity, and can be assumed to remain constant at all red- 
shifts. Based on this result, we built a Hubble diagram 
for quasars, which extends up to z > 6, and is in per¬ 
fect agreement with the analogous Hubble diagram for 
supernovae in the z ~ 0.01 — 1.4 range in common. The 
main advantage of this method is clearly the possibility 
of testing the cosmological model, and measuring the cos¬ 
mological parameters, over a wider redshift range than 
any other cosmological pr obe, with the possible excep- 
tion of gamma-ray bursts (Ghirlanda et al.||2006). The 
main limitation lies in the large observed dispersion of 
the relation (cr^O.3 in a logLuv-logTx plane). As a 
consequence, large samples are needed to obtain signifi¬ 
cant constraints on cosmological parameters. 

With present data, we were able to build a sample of 
808 quasars by merging several literature samples accu¬ 
rately cleaned from BAL, radio loud, and optically heav¬ 
ily reddened objects. A correction for dust extinction was 
applied to the moderately reddened ones. This sample 
allows a first application as a cosmological probe: as¬ 
suming the ACDM model, we obtained HM=0.2lto;?oi 
and a looser constraint on Ha (a lower limit of ~ 0.7, 
due to small amount of low-z quasars in our sample). 
The constraints on the cosmological parameters are ob¬ 
tained by fitting a Hubble Diagram, analogously to what 
is done with supernovae. Therefore, we performed joint 
fits of our data and the Union 2.1 SNe sample of the Su¬ 
pernovae Cosmology Project, and obtained much tighter 


constraints: HM=0.2lto;o4 ^A=0.74^Qgg. These 
uncertainties are still large compared to the estimates ob¬ 
tained by combining all the available cosmological probes 
(CMB, SNe, BAO, leasing). However, they are better by 
about a factor of two that those from SNe alone, show¬ 
ing that the quasar sample available today can already 
provide a significant contribution to cosmological studies 
that are based on distance measurements only. 

Finally, we discussed the potential of quasars as cosmo¬ 
logical probes considering future samples with both UV 
and X-ray measurements, which will become available in 
the next few years. With currently available observato¬ 
ries, new dedicated observations of well-selected high-z 
quasars will greatly improve the test of the cosmolog¬ 
ical model at z > 4. The forthcoming eROSITA all¬ 
sky X-ray survey will provide X-ray measurements for 
more than 100,000 SDSS quasars. Further in the future, 
surveys from Euclid and LSST in the optical-UV, and 
Athena and other possible wide field X-ray survey tele¬ 
scope, will provide samples of millions of quasars. With 
these samples it will be possible to obtain constraints 
on possible deviations from the standard cosmological 
model, which will rival and complement those available 
from the other methods. 

In future papers we will further investigate the poten¬ 
tial of the new method presented here. We will further 
discuss possible limitations and systematic effects, which 
are likely not relevant with relatively small samples such 
as the one available at present, but may become signif¬ 
icant when the increased size of the samples will allow 
more precise measurements. 
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APPENDIX 

APPENDIX A 
LITERATURE SAMPLE 


We have c onsi dered the samples p resented by Steffen et al. (20061, Shemmer et al. (2006), Just et al. (2007), Young 
et al. (2010), and|Lusso et al. (2010) with both optical and X-ray luminosities at 2500 A and 2 keV, respectively. Given 


that uncertainties for both optical and X-ray luminosities were not published in most of these works, we retrieved, 
where possible, all multi wavelength information and re-compute luminosi ties with their uncertainties. For all catalog 
correlations we have used the Virtual Observatory software TOPCAT (Taylor 2005) available onlin^ The total 
quasar sample considered consists on 1,138 objects. A summary of the total quasar sample is given in Table 

A.l. Optically selected samples 


The Steffen et al. samp le (333 sources) contains 155 objects from the SDSS-DR2 quasar catalog (Abazajian et ah, 
2004 Strateva et al.|20(i5 ). For these sources we have updated the optical values using the more recent quasar catalog 


^ http://www.star.bris.ac.uk/~mbt/topcat/ 




































A Hubble diagram for quasars 


13 


published by Krawczyk et al. (2013). We obtained 133 matches (excluding BAL and radio-loud quasars) using a 
matching radius of 3 arcsec. I'hese 133 quasars have been then cross correlated with the the ROSAT archive (104 
matches with a radius of 22 arcsec). The X-ray monochromatic flux at 2 keV has been obtained by converting the 
ROSAT/PSPC count rates in the energy band 0.1-2.4 keV into unabsorbed 0.5-2 keV huxes by using Webpimms 
assuming a power law spectrum with no intrinsic absorption and a photon index P = 2.0 modified by Galactic 
absorption ( Kalberla et al.|[2005 ). 

We adopted the original X-ray data for the remaining sources, while we tried to update the optical lumin osities 
whenever p ossible by cross -matching them w ith the available SDSS catalogs. We considered the SDSS-DR7 (Shen 
et al. 2011), BOSS-DRIO (Paris et al. 2014), and the SDSS catalogs available in TOPCAT. This SDSS sample of 
104 quasars is combined with moderate-luminosity AGNs from the COMBO-17 survey (52 objects), a subsample of 
sources from the Bright Quasar Survey (BQS) quasar catalog (46 objects), low-redshift optically selected AGNs (24 
Syls), and additional optically selected, z > 4 AGNs (54 sources). A summary of the [0.5 — 2] keV limit for each 
subsample is provided by Steffen et collaborators in their Table 3. The new Steffen et al. sample is so composed by 
280 objects. 

We adopted a similar procedure for the quasars presented by|Shemmer et al.| (|2006|) and Just et al. (2007). We 
retrieved all multi wavelength information from infrared (WISE, UKIDSS Data Release 9, and 2MASS) to optical 
(SDSS-DR9) to estimate optical luminosities with their uncertainties, while we kept the original X-ray luminosities 
for both subsamples. 

The sample presented by Just et al. (2007) is originally composed by 34 sources. Their s ample of highly luminous 
quasars has been mainly drawn from the SDSS DR3 quasar catalog (Schneider et al. 2005). Eleven quasars already 
had detections in archival X-ray data, while 21 were observed using the Advanced GGD imaging Spectrometer (ACIS) 
during the Chandra Cycle 7 Guaranteed Time Observing program. The requested snapshot exposure for each target 
was 4 ks, and all targets were strongly detected with ~ 10 — 150 counts from 0.5 to 8 keV (see their Section 2.2 for 
details). Their SDSS sample thus includes 32 quasars with absolute magnitudes values of —29.28 to —30.24, all of 
which have sensitive X-ray coverage; and span the range redshift 1.5 < z < 4.5. Two complementary z > 4 quasars 
were added to their sample APM 08279-1-5255 at z = 3.91 and HS 1603-1-3820 at 2.51, leading to their “core” sample 
of 34 objects. From this core sample we excluded 7 objects classified either BALs, RLQs, and/or lensed sources (see 
their Table 1). We have also neglected the object 142123.98-1-463317.8 at z = 3.36 which has been classified BAL by 
Gibson et al.|(|2009), and the sources 170100.62-1-641209.0 at z = 2.74 and 152156.48-1-520238.4 at z = 2.19 which were 
already in the btetnen et al. sample. The final Just et al. sample considered in our analysis is thus composed by 24 
o bjects. _ 

Shemmer et al. (2006) presented Chandra observations of 21 z > 4 quasars mainly selected from the SDSS DR3 


quasar catalog for z < 5.4 and with near-infrared imaging and spectroscopy for higher redshifts. Nineteen quasars were 
targeted with Chandra during Cycles 4 and 6 with short (3 — 30 ks) X-ray observations. Three other weak (or absent) 
emission-line quasars (WLQs) from DR3 were included in their sample (see their Section 2 for details), but given 
the peculiarity of such objects we neglected them for our study. From their sample we also excluded one radio loud 
quasar (SDSS J001115.23-I- 144601.8 at z = 4.97), one moderate radio emitting quasar (SDSS J144231.72-1-011055.2 
at z = 4.51, also classified WLQ), and two BALs (SDSS J104845.05 -f 463718.3, SDSS J165354.61 -f 405402.1). The 
final Shemmer et al. high-redshift sample considered in our analysis is thus composed by 14 objects. 

We have further increased our sample with the one published by [Young et al. (2010), which is composed by 327 
quasars (their SPECTRA sample) selected by cross-correlating the SDSS DR5 quasar catalog with the XMM-Newton 
archive. They fitted three models to each XMM-Newton spectrum: (1) a single power law with no intrinsic absorption, 
(2) a fixed power law with intrinsic absorption left free to vary, and (3) an intrinsically absorbed power law, with both 
F and JVh left free to vary. Any spectrum without a good fit, or with significant contribution from a strong soft 
excess component or absorption is excluded from the final sample. All sources have both optical and X-ray spectra 
with an X-ray signal-to-noise ratio (S/N) > 6 spanning a redshift of z = 0.1 — 4.4 (with an j—band magnitude of 
15.2 — 20.4, see their Section 2 for further details). We have neglected overlapping objects leading to a final sample 
of 278 sources. Optical luminosities have been updated for 242 objects by cross-matching them with the SDSS-DR7 
quasar catalog using a matching radius of 3 arcsec. For the remaining sources (36 quasars) we considered the published 
optical luminosities. We also kept the X-ray luminosity values as published by Young et collaborators. 

A.2. X-ray selected sample 


We considered an updated version of the catalo g already published by Brusa et al. ( 2010[ ), which includes the 
photometric redshift catalog by Salvato et al. (2011), and new spectroscopic redshift measurement^ We have selected 
1375 X-ray sources detected in the 0.5-2 keV band at a flux larger th an 5x 10~^^erg s~^ cm~^ over the COSMOS area (2 
deg^) and for which a reliable optical counterpart can be associated (Brusa et al. |2010[ ). From this sample, 426 objects 
are spectroscopically classified as broad-line AGN on the basis of broad emission lines {FWHM > 2000km s“^, 
“spectro-z” sample hereafter) in their optical spectra. In order to extend our XMM-COSMOS sample to fainter 
magnitudes, we added to the spectro-z sample a sample of 116 Type-1 AGN defined as such via S ED-fitting (“photo- 
z” sample hereafter). The photo-z sample has been selected following the same approach as in Lusso et al. (2013 


® The multi-wavelength XMM-COSMOS catalog can be re- 
trieved from: 

http://www.mpe.mpg.de/XMMCosmos/xmm53-release/, version 


fst November 2011. 
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TABLE 4 

Summary of the total quasar sample. 


Number [0.5 — 2] keV limit Area z 

erg sec“^ cm“^ deg^ 


280 

24 

14 

542 5 X 10-1® 

278 


0.009- 6.280 
1.760- 4.610 
4.720- 6.220 
2 0.041 - 4.255 

0.160 - 4.441 


1,138 0.009-6.280 


i mag Reference 


12.20 - 23.70 

Steffen et al. ( 

2006 1 

15.00 - 20.20 

Just et al. 

20071"^ 

18.34 - 23.78 

Shemmer et al. 

12^0061 

16.86 - 26.04 

Lusso et al. i 2U10ll 


15.26 - 20.40 

Young et al. (|2010 



15.00 - 23.78 Total 


‘‘Not well defined. 


see their Section 2 for details). We selected all sources with a be st-fit photometric cl assification consistent with an 
AGN-dominated SED (i.e., 19 < SED — Type < 30 as presented by|Salvato et al.||2009|). The X-ray luminosity values 
have been computed following the same approach as outlined in Section 2.1 in Lusso et al. (2013). Briefly, count rates 
in the 0.5-2 keV and 2-10 keV are converted into monochromatic X-ray fluxes in the observed frame at the geometric 
mea n of the soft (1 keV) a nd hard (4.5 keV) energy bands using a Galactic column density A^h = 2.5 x 10^° cm“^ 
(see Gappelluti et al.|2009 1, and assuming a photon index T = 2 and T = 1.7, for the soft and hard band, respectively. 
Fluxes are then blueshitted to the rest-frame. The rest-frame monochromatic flux at 2 keV is finally obtained by 
interpolation of these fluxes if the source redshift is lower than ^1, by extrapolation considering the slope between the 
fluxes described above for higher redshifts. The XMM-GOSMOS sample comprises 542 AGN (426 with spectro-z and 
116 with photo-z) spanning a redshift range of z = 0.041 — 4.255. 

A.3. Rest-frame optical luminosity 

To obtain the rest-frame monochromatic luminosities at 2500A we used all the available photometry compiled in 
the SDSS and XMM-GOSMOS catalogs. To compute the rest-frame AGN SEDs we considered the flux density 
(erg/cm^/s/A) and 1 -ct uncertainty at the effective wavelength of the filter listed in the catalogs. Galacti c reddening 
has be en taken into account: we used the selective attenuati on of the stella r conti nuum k{X) taken from Fitzpatrick 
(1999) with Ry = 3.1. Galactic extinction is estimated from Schlegel et al. (1998) for each object. We derived total 
luminosities at the rest-frame frequency of the filter according to the standard formula 

(Al) 


= v„Fi, AttD 


L- 


The data for the SED computation from mid-infrared to UV (upper limits are not considered) were then blueshifted 
to the rest-frame and no K-correction has been applied. We determine a “first order” SED by using a first order 
polynomial function, which allows us to build densely sampled SEDs at all frequencies. This choice is motivated by 
the fact that a single interpolation with a high-order polynomial function could introduce spurious features in the 
final SED. In the case 2500 A is covered by no less than two data points, the Tuv values are extracted from the 
rest-frame SEDs in the logv — \ogvL^ plane {vFy oc iF). If the SED is constructed by two data points (or more), but 
they do not cover 2500 A, luminosities are extrapolated by considering the last (first) two photometric data points. 
Finally, we corrected the Ruv estimates taking into account the redshift-dependent contribution of emission lines to 
the photometric points, as described in Appendix A.4. Uncertainties on monochromatic luminosities (L^ oc v~^) from 
interpolation (extrapolation) between two values Li and L 2 are computed as 




(<5L2)A 


(A2) 


A.f- Emission lines contamination 

Broad-band photometry of quasars usually contains a certain line contribution. To quantify th e amplitude of this 
contribution we have compared the continuum flux densi ty at rest-frame 2 500 A as compile d by 


with our estimates obtained as described in appendix A.3 Optical fluxes in Shen et al. (2011) were obtained through 


Shen et al. 


( 2011 ) 


a fit of the SDSS spectra where five parameters were simultaneously fitted: fhe normalization and slope of the power- 
law continuum, and the normalization, line broadening and velocity offset relative to the systemic redshift for the 
iron template fit. The interested reader should refer to their Section 3 for details. We note that continuum flux 
measurements in the Shen et al. catalogue were neither corrected for intrinsic extinction/reddening, nor for host 
contamination, while our UV fluxes are corrected for dust extinction only. How ever, our sample selection reduces both 
reddening and host contaminations at minimum as discussed in Section |3.1| The left panel of Figure shows the 
comparison of the two Fuv measures (where both are available) for 448 objects within our selected sample. 

Our flux estimates are in good agreement with the ones c omputed by Shen et ah, yet systematically higher. The 
average values of the difference between the optical flux by Shen et al. (2011) and our measurements (AlogEuv = 
logFuv,Shen+ii — logEyv) are plotted as a function of the average redshift iii each bin. The solid and dashed lines 
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- 10 ][ergs ' cm 2] 


_ Fig. 9.— Left panel: comparison between the optical flux measured using described in appendix |A.3| with the flux measured by |Shen| 

|et al.| (|2011|l through a complete spectral fit. Colors and point shapes represent different redshift intervals, as indicated in the plot. I'he 
dastied red line re presents the one-to -one relation. The inset plot shows the average values (along with their errors) of the difference between 
the optical flux by|Shen et aL|(|2011|l and our measurements (A log ihjv = log ^uv,Shen+li ~ logTuv) as a function of the average redshift 
in each bin (chosen to have approximately the same number of sources). The solid and dashed lines represent the mean and the error on 
the mean of the unbinned A log -Puv distribution (mean is ~ —0.05 with a = 0.12). Right panel: same for X-ray fluxes in the 2-10 keV 
energy band. The comparison is between our estimates and those from the spectral analysis in Mainieri et al. (2011). The mean over the 
whole sample is -0.02 with a = 0.29. 


represent the mean and the error on the mean of the unbinned Alogidjv distribution (mean is ^ —0.05 with cr = 0.12 
in logarithm). A correction to the values of idjv has been applied to our data, with no significant change in the results. 
This effect has been also considered in the simulations, where we found that with large samples allowing more precise 
measurements of the cosmological parameters, the effect is not negligible (see Section 7.2). 


A.5. X-ray K—correction 

For 646 quasars (~57% of the total sample) we have adopted an X-ray AT—correction, whose systematics may affect 
the X-ray flux measurements and, in turn, our cosmological results. To quantify the amplitude of this possible effect, 
we have considered a subsample of quasars within the XMM-COSMOS sample for which the rest-frame 2-10 keV fluxes 
were available from a spectroscopic analysis. Thes e flu xes are then compared with the ones computed by integrating 
t he rest-frame X-ray S EP as estimated in Section [A.2| in the 2-10 keV range. 

Mainieri et al.| (|2011 1 presented a detailed spectral analysis of X-ray sources in the XMM-Newton COSMOS field. 


Here we consider a sample of 408 quasars with more than 70 net counts in the 0.3-10 keV energy band. The rest frame 
2-10 keV fluxes (F'[ 2 _io]keV.sp) are estimated from a fit of the XMM-Newton spectra using an intrinsically absorbed 
power law with both T and (Vh free to vary. The power law is then extrapolated to lower energies considering the best 
fit photon index, and the Arg^iojkeV.sp is finally estimated correcting the values for intrinsic absorption. 

The ri ght panel of Figure s hows the comparison between the rest-frame 2-10 keV fluxes from the spectral analysis 
done by Mainieri et al. (20111 as a function of our 2-10 keV flux measurements. We considered this correction in 
both our data and in the simulations. In the present sample, the effect is totally negligible. Regarding the simulated 
samples, the deviations found here may in principle alter the measurement of the cosmological parameters. However, 
the future X-ray surveys will allow a direct spectral analysis for all sources, so such systematics in the estimates of Fx 
will be absent. 


APPENDIX B 
TESTS OE THE METHOD 

We tested the method used to derive the cosmological parameters, based on the Hubble diagram of quasars, with 
simulated quasar samples and different values of and Ha. We started from the Fuv ~ -2 distribution of our sample 
of 808 quasars, and performed the following simulations: (1) we simulated Fx from (z,Fuv) of each quasar in our 
sample, assuming four different combinations of Hm and Ha: (Hm,Ha)=( 0.2,0.2), (0.2,0.8), (0.8,0.8), (0.8,0.2); (2) we 
repeated the simulations for samples with the same (z,Fuv) distribution, but made of 50,000 objects. In this way, we 
test the ability of measuring different parameters of our actual sample, and also the correctness of the method in the 
limit of very large samples. In all cases we assume a dispersion of the logFyv-logFx relation <5=0.30. The results are 
shown in Figure 9, and demonstrate the solidity of our approach. 

REFERENCES 

Abazajian, K., Adelman-McCarthy, J. K., Agueros, M. A., et al. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 

2004, AJ, 128, 502 219, 12 

Aird, J., Comastri, A., Brusa, M., et al. 2013, ArXiv e-prints Avni, Y., & Tananbaum, H. 1986, ApJ, 305, 83 










































16 


Risaliti & Lusso 



Qm 

Fig. 10.— Contours obtained fitting the Hubble diagram of quasars for simulated samples with the same size as the “real” one (empty 
contours, 68% confidence level), and for samples of 50,000 objects (filled contours, 68% and 95% confidence levels), and assuming four 
different pairs of values for flM and Ha, as indicated by the black stars. 


Baldwin, J. A. 1977, ApJ, 214, 679 

Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, A22 

Brandt, W. N., Laor, A., & Wills, B. J. 2000, ApJ, 528, 637 

Brusa, M., et al. 2010, ApJ, 716, 348 

Cappelluti, N., et al. 2009, A&A, 497, 635 

Edelson, R., & Malkan, M. 2012, ApJ, 751, 52 

Edelson, R., Gelbord, J. M., Horne, K., et al. 2015, ArXiv e-prints 

Fitzpatrick, E. L. 1999, PASP, 111, 63 

Foreman-Mackey, D., Hogg, D. W., Lang, D., Sz Goodman, J. 

2013, PASP, 125, 306 

Gallagher, S. C., Brandt, W. N., Sambruna, R. M., Mathur, S., & 
Yamasaki, N. 1999, ApJ, 519, 549 
Ghirlanda, G., Ghisellini, G., & Firmani, C. 2006, New Journal of 
Physics, 8, 123 

Gibson, R. R., Brandt, W. N., Sz Schneider, D. P. 2008, ApJ, 685, 
773 

Gibson, R. R., Jiang, L., Brandt, W. N., et al. 2009, ApJ, 692, 758 
Gilli, R., Vignali, C., Rosati, P., et al. 2010, in COSPAR Meeting, 
Vol. 38, 38th COSPAR Scientific Assembly, 2361 
Green, P. J., Schartel, N., Anderson, S. F., et al. 1995, ApJ, 450, 
51 

Grupe, D., Komossa, S., Leighly, K. M., & Page, K. L. 2010, 
ApJS, 187, 64 

Hao, H., Elvis, M., Givano, F., et al. 2014, MNRAS, 438, 1288 
Hopkins, P. F., Strauss, M. A., Hall, P. B., et al. 2004, AJ, 128, 
1112 

Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90 
Ivezic, Z., Tyson, J. A., Abel, B., et al. 2008, ArXiv e-prints 
Just, D. W., Brandt, W. N., Shemmer, O., et al. 2007, ApJ, 665, 
1004 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, 
A&A, 440, 775 

Kilerci Eser, E., Vestergaard, M., Peterson, B. M., Denney, K. D., 
Sz Bentz, M. C. 2015, ApJ, 801, 8 
Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 
330 

Krawczyk, C. M., Richards, G. T., Mehta, S. S., et al. 2013, 
VizieR Online Data Catalog, 220, 60004 
La Franca, F., Bianchi, S., Ponti, G., Branchini, E., Matt, G. 

2014, ApJ, 787, L12 


Laureijs, R., Gondoin, R, Duvet, L., et al. 2012, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 8442, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, 0 
Lusso, E., et al. 2010, A&A, 512, A34 

Lusso, E., Hennawi, J. F., Comastri, A., et al. 2013, ApJ, 777, 86 
Mainieri, V., Bongiorno, A., Merloni, A., et al. 2011, A&A, 535, 
A80 

Melia, F. 2014, J. Cosmology Astropart. Phys., 1, 27 
Merloni, A., Predehl, P., Becker, W., et al. 2012, ArXiv e-prints 
Mortlock, D. J., Patel, M., Warren, S. J., et al. 2012, MNRAS, 
419, 390 

Murray, S., Gilli, R., Tozzi, P., et al. 2009, in Astronomy, Vol. 
2010, astro2010: The Astronomy and Astrophysics Decadal 
Survey, 217 

Osmer, P. S., Sz Shields, J. C. 1999, in Astronomical Society of 
the Pacific Conference Series, Vol. 162, Quasars and 
Cosmology, ed. G. Ferland Sz J. Baldwin, 235 
Page, M. J., Brindle, G., Talavera, A., et al. 2012, MNRAS, 426, 
903 

Paris, I., Petitjean, P., Aubourg, E., et al. 2014, A&;A, 563, A54 
Phillips, M. M. 1993, ApJ, 413, L105 

Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2015, 
ArXiv e-prints 

Prevot, M. L., Lequeux, J., Prevot, L., Maurice, E., Sz 
Rocca-Volmerange, B. 1984, A&A, 132, 389 
Richards, G. T., et al. 2006, AJ, 131, 2766 

Risaliti, G., Salvati, M., Sz Marconi, A. 2011, MNRAS, 411, 2223 
Salvato, M., et al. 2009, ApJ, 690, 1250 
—. 2011, ApJ, 742, 61 

Schlegel, D. J., Finkbeiner, D. P., Sz Davis, M. 1998, ApJ, 500, 
525 

Schneider, D. P., Hall, P. B., Richards, G. T., et al. 2005, AJ, 

130, 367 

Shemmer, O., Brandt, W. N., Schneider, D. P., et al. 2006, ApJ, 
644, 86 

Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 
45 

Steffen, A. T., Strateva, L, Brandt, W. N., et al. 2006, AJ, 131, 
2826 

Strateva, 1. V., Brandt, W. N., Schneider, D. P., Vanden Berk, 

D. G., & Vignali, C. 2005, AJ, 130, 387 








A Hubble diagram for quasars 


17 


Suzuki, N., Rubin, D., Lidman, C., et al. 2012, ApJ, 746, 85 
Tananbaum, H., Avni, Y., Branduardi, G., et al. 1979, ApJ, 234, 
L9 

Taylor, M. B. 2005, in Astronomical Society of the Pacific 
Conference Series, Vol. 347, Astronomical Data Analysis 
Software and Systems XIV, ed. P. Shopbell, M. Britton, & 

R. Ebert, 29 

Vagnetti, F., Antonucci, M., V Trevese, D. 2013, A&A, 550, A71 
Vagnetti, F., Turriziani, S., Trevese, D., &; Antonucci, M. 2010, 
A&A, 519, A17 


Wang, J.-M., Du, P., Hu, C., et al. 2014, ApJ, 793, 108 
Watson, D., Denney, K. D., Vestergaard, M., V Davis, T. M. 
2011, ApJ, 740, L49 

Young, M., Elvis, M., & Risaliti, G. 2010, ApJ, 708, 1388 
Zamorani, G., Henry, J. P., Maccacaro, T., et al. 1981, ApJ, 245, 
357 



